Optical conductivity in the CuO double chains of PrBa 2 Cu 4 08: 
Consequences of charge fluctuation 
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We calculate the optical conductivity of the CuO double chains of PrBa2Cu4 0s by the mean- 
field approximation for the coupled two-chain Hubbard model around quarter filling. We show that 
the ~40 meV peak structure, spectral shape, and small Drude weight observed in experiment are 
reproduced well by the present calculation provided that the stripe-type charge ordering presents. 
We argue that the observed anomalous optical response may be due to the presence of stripe- 
type fluctuations of charge carriers in the CuO double chains; the fast time scale of the optical 
measurement should enable one to detect slowly fluctuating order parameters as virtually a long- 
range order. 

PACS numbers: 71.10.Fd, 73.90.+f, 71.27.+a, 78.20.Bh, 71.30.+h 
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I. INTRODUCTION 

The CuO double chains of PrBa 2 Cu40 8 (see Fig. 1(a)) 
show metallic conductivity along the chain direction 
down to 2 K Jl], . Because the system shows insulating 
conductivity along the directions perpendicular to the 
chains (except at low temperatures <140 K), and also 
because the holes in the Cu02 planes are localized in the 
O 2p-orbitals around Pr ions §, PrBa 2 Cu 4 8 may be re- 
garded as a model material of quasi-one-dimensional (ID) 
correlated conductors except at low temperatures. A 
possible Tomonaga-Luttinger-liquid description has then 
been tried to explain available experimental data ^] . 

A number of anomalous behaviors have however been 
reported in the charge degrees of freedom of this sys- 
tem a nuclear-quadrupole-resonance (NQR) ex- 
periment has demonstrated the presence of anomalous 
spin-lattice (1/Ti) and spin-spin (1/T 2 ) relaxations || 
and optical conductivity experiment clearly shows the 
presence of anomalous charge excitations at w > 40 meV 
pj. The CuO double-chain system therefore provides a 
good opportunity for studying the anomalous charge dy- 
namics of strongly correlated ID electron systems, such 
as charge ordering (CO) and charge fluctuations, just as 
the quasi- ID organic systems do Q, For example, a 
possible mechanism of metallization of this system due 
to charge frustration has recently been proposed 0. 

Motivated by such development in the field, we in this 
paper focus on the optical conductivity spectra cr(u>) ob- 
served by Takenaka et al. [Q and consider the origin of the 
appearance of the characteristic peak structure at w ~ 40 
meV as well as the presence of extremely small Drude 
weight. We use the mean- field approximation to the ex- 
tended Hubbard model defined on a lattice of coupled 
two chains to obtain the ground-state phase diagram, 
where we find that there appear two types of CO phases, 
depending on the parameter values. We then calculate 
the optical conductivity spectrum a{u>) for each of these 
two phases in the mean-field approximation. We find 
that the spectra calculated for the ground state of the 



stripe-type charge ordering reproduces the experimental 
peak structure including not only its peak position and 
spectral shape but also the observed very small Drude 
weight. We also find that the spectra calculated for the 
in-line-type CO phase completely fails to reproduce the 
experimental spectra. 

It is known that the CuO double chains do not show 
the long-range CO but have the fluctuation of electronic 
charges with very slow time scale of the order of 0.1 — 10 
MHz ||. We then point out that the fast time scale 
of the optical measurement should enable one to detect 
slowly fluctuating order parameters as virtually a long- 
range order; i.e., the spectra observed should basically be 
the same as those for the long-range CO state. There- 
fore, we argue that the fact that our mean-field optical 
spectra agree well with experiment reasonably suggest 
that the stripe-type charge fluctuations should present 
in the system and should be the origin of the ^40 meV 




FIG. 1: Schematic representations of (a) the CuO double 
chain and (b) corresponding two-chain Hubbard model. Note 
that the :r-axis (y-axis) is taken along the crystallographic 
6-axis (a-axis). The parameters t and V (t' and V) are as- 
sociated with the chain (zigzag) bonds. The stripe-type CO 
discussed in the main text is schematically shown as shaded 
circles in (a). 
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peak structure in the optical conductivity spectrum of 
the CuO double chains of PrBa 2 Cu40 8 . 

This paper is organized as follows. In Sec. II, we define 
a model for describing the electronic states of the CuO 
double chains. In Sec. Ill, we derive the mean-field phase 
diagram of the Hamiltonian, and in Sec. IV, we calculate 
the optical conductivity of the model. The results are 
discussed in Sec. V by comparing with experimental data 
and their implication is considered. Summary is given in 
Sec. VI. 



II. MODEL 

The structure of the CuO double-chain system is illus- 
trated in Fig. 1(a), where we show the two CuO chains 
connected with zigzag bonds. The electronic structure 
may be modeled as a single-band extended Hubbard- 
type Hamiltonian representing the doping holes in the 
antibonding band of the Cu 3d x 2_ y 2 and O 2p a or- 
bitals. The filling of holes in the chains is reported to be 
ri ~ 0.5 (quarter filling) [|, ||; the angle-resolved photoe- 
mission spectroscopy experiment Q suggests the value 
n ~ 0.46 ± 0.06 and optical conductivity measurement 
H suggests the value n ~ 0.4. 

We therefore adopt in this paper the single-band Hub- 
bard model around quarter filling defined on the lattice 
of two chains connected with zigzag bonds as shown in 
Fig. 1(b). The Hamiltonian may be of the following form: 

H= - *tf( c L<> +H.c.) + U^2nifnu 

<ij><J i 
+ VijUitlj (1) 

<y'> 

where cj^ (ci a ) is the creation (annihilation) operator of 
a hole at site i and spin a (=|, J,), and rii ~ m\ + 
is the number operator. We take into account the hop- 
ping parameters for the chain bond t (— i 2 ) and for the 
zigzag bond t 1 (= ti) where we note t' <t because of the 
bonding angle minimizing the overlap of the Cu 3d x 2_ y 2 
and O 2p a wave functions for the hopping parameter t'; 
in the following, we use the value t'/t = 0.1 confirm- 
ing that the sign of t' does not change the results. Wc 
also take into account the intersite Coulomb interactions 
V for the chain bonds and V for the zigzag bonds, as 
well as the on-site Coulomb interaction U; the values 
are varied to simulate various situation. We restrict our- 
selves to the case at and less than quarter-filling of holes 
n = N/L < 0.5 where N and L are the total numbers of 
holes and lattice sites, respectively. 

III. MEAN-FIELD PHASE DIAGRAM 

We adopt the mean-field approximation 

n lT n u ~ (n^)nn +n iT (n u ) - (n lT )(n a ) (2) 



for the interaction terms of the Hamiltonian Eq. (1). We 
assume periodicity of the ordered states where each unit 
cell r contains M sites; the lattice site may thus be spec- 
ified as i — (r, /i) with the unit cell r (= 0, • • • , L/M — 1) 
and site position /j, (= 1, • • • , M) therein. We have 2M 
order parameters 



s» = (rvr) + Ku) - n 
2^v = M - M 



with the charge neutrality constraint 5^ = 
define c u ka as the Fourier transform of Ci„ = c r 



L/M 

^ e ikRr »c 



rfia 



(4) 
(5) 

0. We 



(6) 



where R r ^ = rMa + p,a and k is the wavevector taking 
L/M values in the Brillouin zone —ir/Ma < k < ir/Ma. 
a is the lattice constant of the topologically equiva- 
lent single-chain model with nearest-neighbor (t' , V') and 
next-nearest-neighbor (t, V) interactions. The mean-field 
Hamiltonian is then written as 



MF 



H 



with 



where 



E 

ho 



(c f c< 



'MkaJ 



( Cik 

C2k 



\ CMka J 



[Hfj\ = [H t ] + [H}*\ 



(7) 



(8) 



if v = (j, ± 1 
if v = fi ± 2 
if n= = M 
if H = l,v = M 
if ^ = 2, i/ = M 
if n = M -1,1/ 
if n = M, v = 1 
if n = M, v = 2 
otherwise 



- 1 



= 1 



(9) 



and 



H 



MFl 



U{n„,- a ) + ^((n^-i) + (n„+i)) 



+V((n^ 2 ) + (n M+2 )) 



(10) 



In Eq. (9), we define a x 
have 



-U e 



±ilka 



(1 = 1,2). Also, we 



E r = - 



M 



Uifij ~ (n,-)rij + ni(rij) - (ni)(rij) 



(3) 



(V'in^in^a) + V{7i^)(n^ +2 ,a)) (-11) 
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Diagonalizing the matrix (8) by the unitary transforma- 
tion 

[^"(pS.pS.-.P^) (12) 

•4.1. ( (")* («)* M *\ , 

with p k J = (p\ k > a ,p\ k > a , ■ • ■ ,pVkJ> we have 

^ MP = 5:^7^7^+^ (13) 

nfccr 

where E k J is the energy of the n-th quasiparticle. The 
annihilation operator of the quasiparticle may be written 
as 

M 

Ink* = J2 p( ?k* C ^ ( W ) 
/i=l 

M 

C^ka = ^P^Llnka- (15) 

71=1 

The mean-field ground state is then written as 

occ. 

i°) = n T«fc CT i vacuum >- (i 6 ) 

nka 

The order parameters are determined self-consistently so 
as to minimize the total energy of the system. The rela- 
tion 

<^) = ^EI*r^-4: ) ) (17) 

nk 

is used to make iterations for self-consistency where 9{x) 
is the Fermi function at T = K and £p is the Fermi 
level. 

The calculated ground-state phase diagram of the 
model is shown in Fig. 2, where we use M = 8. We find 
that there appear the following three types of phases: (i) 
The charge-ordered phase denoted as CO-1 where the in- 
line-type ordering of charges emerges (see inset of Fig. 2). 
The spins are aligned antiferromagnetically at sites on 
one of the two chains; there is no spin polarization on the 
other chain (and even charges vanish if we set t' = 0). 
(ii) The charge-ordered phase denoted as CO-2 where 
the stripe-type ordering of charges emerges (see inset of 
Fig. 2). The up and down spins on the neighboring two 
sites would form a spin-singlet pair if quantum fluctu- 
ation is introduced, (iii) The spin-density-wave (SDW) 
phase with a small modulation of the charge-density dis- 
tribution, which we denote as SDW+CDW, as shown in 
Fig. 2. 

We first note that the U/t dependence of the phase 
boundary is rather small, that the CO phases are main- 
tained up to ~10% doping levels, and that the doping 
dependence of the location of the phase boundary is also 
very small. Because we assume the unit cell of 8 sites 
with the commensurate CO, the system becomes metal- 
lic upon doping where the small Fermi surface emerges. 
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FIG. 2: Mean-field phase diagram of the two-chain Hubbard 
model at quarter filling (n = 0.5, solid line with filled circles) 
and at a 10% doping level (n = 0.45, dashed line with open 
circles). Two types of the charge ordering CO-1 (in-line-type) 
and CO-2 (stripe-type) are illustrated schematically in the 
figure, where the circles with an arrow indicate the amount 
of charge carriers localized on the sites and arrows indicate 
the direction of the spins of carriers. When t' is smaller than 
a critical value, the pure SDW phase appears around V = 
V' = 0. This region changes into the SDW+CDW (weak 
CO-2) phase when t' becomes larger than the critical value 
(as shown in the figure) . This phase is continuous to the CO-2 
phase; the hatched region indicates the rapid crossover from 
the SDW+CDW phase (for small V) to the CO-2 phase (for 
large V). 



We also note that the phase boundary between CO-1 
and CO-2 is discontinuous; i.e., the transition is of the 
first order. The transition between the phases CO-2 and 
SDW+CDW is however continuous; when t' is smaller 
than a small critical value, a pure SDW phase is realized 
in this 'SDW+CDW region and the transition to the 
CO-2 phase is of the second order. When t' exceeds the 
critical value, the phase boundary becomes obscure and 
the transition changes into a crossover; when the value 
of U/t decreases, the crossover becomes broader. 

The spin degrees of freedom of the fluctuating CO 
phases seems intriguing; the CO-1 phase, if it is a long- 
range CO, should behave like a ID Heisenberg antifer- 
romagnet, and in the CO-2 phase, the spin-singlet for- 
mation should occur. Experimentally, it has been found 
|S that the temperature- independent Pauli paramagnetic 
contribution and Curie- Weiss contribution coexist with a 
small anomaly due to the ordering of Pr spins at low tem- 
peratures. It would be interesting if any influence of the 
fluctuating charges on the spin states can be detected in 
experiment. 
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IV. OPTICAL CONDUCTIVITY 



with 



We now calculate the optical conductivity spectra for 
the states derived in the previous section. Using the 
eigenstates |/) and energies Ef (where / = denotes 
the ground state) of the many-body Hamiltonian Eq. (1), 
we generally define the optical conductivity a aa {oj) (a = 
x,y,z) as 

ffaa(w) = D aa 8{u>) 



ire 



^^|(/|ja|0)|V-(S/-^0)I18) 



fro 

with the current operator 



(19) 



where D aa is the Drude weight. In 2D, the sum rule 

MM 

J du>(a xx (u)+a yy (u)) = —(-T) (20) 

holds with the kinetic energy operator 

T = -^t i j\Ti-T j \ 2 clc ja . (21) 



The relation (20) may be used to calculate the Drude 
weight D as 



D 



1 

r 



<\T\0) 



+ 2 y m.\o)\ 2 + 2 y \(f\ jy \o)\\ (22) 

Note that because the present system is finite along the y- 
direction the Drude contribution vanishes for the electric 
field perpendicular to the chains; D is thus associated 
only with E \\ x. 

We use the mean-field results to calculate the optical 
conductivity and Drude weight. The interband optical 
transition from the mean-field ground state may be writ- 
ten 



\S <- t)ka = 7]fc CT 7tfca|0) 



(23) 



where the momentum is conserved. We may then rewrite 
the terms in Eq. (22) as 



E 

fro 

£■ 

fro 



\(f\3x\0)\ 

Ef — Eo 



EE' 

S^t k(7 



*li»|o> 



7 {s) 
J ka 



Ef — Eq ^ E&-eP ~ 



Us - t\m = msin0 J s ( S <?£>(! - *«) (25) 
where 



J. 



M 

\ " / ^ (s)* 

fl+l,ka ^fiko 

H=l 
M 



(Ix) _ ST^ (jp (a)* (t) _ p * (a)* (i) S 



j(iv) _V( 



(-1)" (*i P $?XM p% - Fi pZ ptim) 



(0 

v 



AH 



and 



Ft 



= 0(e F 
It, e lkt 



E 



1,2). 



(27) 



Here, ^ is the angle between the chain and zigzag bonds 
defined in Fig. 1. Also, the kinetic energy term may be 
written 



|T|0> = -a 5 



nkcr 



T 



(i) 



2 cos^ T 



(2) \ g (n) 



(28) 



with 



M 



(0 

nka 



fip ( S:p ( :Im)- (m) 



Note that in the paramagnetic state the spectral weight 
at u> > completely vanishes, leaving only the Drude 
weight, because the present system has only a single band 
and therefore has no interband transitions. 

Calculated results for the parameter and doping de- 
pendences of the optical conductivity spectra at u > 
for E || x are shown in Fig. 3. We note that the spectral 
weight for E \\ y is two orders of magnitude smaller than 
that for E \\ x (and therefore not shown here) because the 
spectra for E || x have the weight roughly proportional 
to t 2 while the spectra for E \\ y have the weight pro- 
portional to t' 2 , and here we assume the value t' ft = 0.1. 
This is consistent with the fact that the experimental 
spectral weight for E \\ a is more than two orders of mag- 
nitude smaller than the weight for E \\ b Q. Conversely, 
this agreement with experiment demonstrates that the 
value of t' is indeed much smaller than the value of t as 
we have assumed in this paper. 

We find in Fig. 3 that the spectra for the CO-2 ground 
state locate in the energy around uj/t — 0.2 — 2 whereas 
the spectra for the CO-1 ground state locate in a much 
higher energy range cj/t = 3 — 6. Because the gaps of the 
quasiparticle band structure of the two CO states have a 
similar value, this difference in the location of the spec- 
tral weight is purely the matrix-clement effects, i.e., the 
selection rule for these two types of CO states is quite 
different. We may have the following real-space inter- 
pretation for the obtained spectra: In the CO-1 ground 
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FIG. 3: Optical conductivity spectra a(ui) with E || x cal- 
culated for the ground states of CO-2 (in (a) and (b)) and 
CO-1 (in (c) and (d)). Drude peak (u) = 0) is not included. 
The results are at quarter filling (n = 0.5) in (a) and (c), and 
at a 2% doping level (n = 0.49) in (b) and (d). Parameter 
values used for calculation are: V'/t = 0.707, V/t = 0.5, and 
t'/t = 0.1 in (a) and (b), and V/t = 1 and t'/t = 0.1 in (c) 
and (d). 



state, the electric field along the chain takes out a hole 
and put it on the neighboring site to make either a doubly 
occupied site via t or a singly occupied site via t', the for- 
mer of which costs the energy U and the latter of which 
costs the energy V'. Since we assume f 3> t' ', the main 
spectral weight appears at the position uj ~ U . We find 
in Figs. 3 (c) and (d) that the position of the main peak 
indeed depends strongly on U but not on V . In the CO-2 
ground state, the excitation by the electric field along the 
chain does not create doubly occupied sites in the strong 
coupling limit, the energy cost of which is at most V. 
We thus have the main peak much lower in energy. We 
find some fine structures of the spectra which reflects the 
fairly complicated quasiparticle multi-band structure. 

Integrated spectral weight of the Drude (u> — 0) and 
interband (uj > 0) contributions are shown separately in 
Fig. 4. We find that the Drude weight grows very rapidly 
by doping away from quarter filling at which the Drude 
weight vanishes. We also find that the U/t dependence 
of the integrated spectral weight is not very strong. 




FIG. 4: Filling dependence (left panel) and U/t dependence 
(right panel) of the spectral weight at u> — (Drude weight, 
filled circles) and integrated weight over uj > excluding the 
Drude contribution (open circles). The CO-2 phase is as- 
sumed. Parameter values used for calculation are: V/t — 0.5, 
V'/t = 0.707, and t'/t = 0.1. 



V. DISCUSSION 

Let us first summarize the experimental aspects of the 
optical conductivity of PrBa2Cu40s observed at room 
temperature 0. The spectra have the following three 
features; (i) a broad peak structure with large spec- 
tral weight (98% of the total weight) located around 
uj ~40 meV, (ii) very small Drude weight (2% of the total 
weight), and (iii) the spectra for E \\ a have the weight 
which is two orders of magnitude smaller the weight for 
E || b. Another important aspect of the charge dynamics 
of this material is given by the NQR experiment [gj ; (iv) 
the CuO double chains have a very slow fluctuation of 
the electric-field gradient caused by the spatial fluctua- 
tion of electronic charge carriers. The system may thus 
be in the vicinity of the long-range CO state. 

Then, in Fig. 5, we compare our calculated opti- 
cal conductivity with the experimental spectrum [|J for 
PrBa2Cu40s. We find that the experimental spectrum 
is well reproduced; (i) The calculated spectrum, which 
is in reasonable agreement with experiment, is for the 
CO-2 ground state near the crossover region to the 
SDW+CDW phase (see Fig. 2). The spectra calculated 
for the CO-1 states completely fail to agree with experi- 
ment as already pointed out in Sec. IV. (ii) The spectra 
for the electric field perpendicular to the chains have the 
weight two orders of magnitude smaller than those for 
the electric field parallel to the chains, which is consis- 
tent with experiment (see Sec. IV). (iii) The calculated 
peak structure is found to locate at rather higher energies 
than experiment. Although this is partly due to insuf- 
ficiency of the mean-field approximation, fine-tuning of 
the parameter values would improve the agreement with 
experiment, (iv) The observed very small Drude weight 
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FIG. 5: Calculated optical conductivity spectrum compared 
with the experimental data of Takenaka et al. [0] (inset). Pa- 
rameter values used for calculation are: t = 0.3 eV, U/t=2.5, 
V/t = 0.5, V'/t = 0.707, t'/t = 0.1, and n = 0.495. The 
Drude peak and the peak at u > have the weight of 6.7% 
and 93.3% of the total weight, respectively. The 5-function of 
the Drude peak is broadened artificially with the Lorentzian 
of width r)/t = 1/IOtt. 



is well reproduced by our calculation, (v) The asym- 
metric spectral shape of the main peak structure, i.e., 
sharp in the lower energy side and broad in the higher 
energy side, is well reproduced by our calculation. The 
width of the peak is however rather narrow compared 
with experiment. This may be due to the effect of elec- 
tron correlations neglected in the mean-held approxima- 
tion, i.e., the scattering of the quasiparticles may well 
broaden the spectra 0] . (vi) As for the amplitude of the 
ordered charges, we find the calculated values of the or- 
der parameters to be 5 = 0.013 (per site) and 2S Z = 0.21 
(per site) for the parameter values used in the calcula- 
tion of the optical conductivity of Fig. 5. This value of 
S is consistent with the experimental estimation where it 
has been pointed out that less than 0.02 holes per Cu are 
frozen 

Finally, let us discuss the optical conductivity spectra 
for fluctuating charge carriers. We should first note that 
there is no experimental indication of the long-range CO 
in PrBa2Cu40g but rather there is only the slow fluctu- 
ation of electronic charge degrees of freedom Then, 
in this situation, we find that the spectra calculated for 
the ordered state agree reasonably well with experimen- 
tal spectra. The time-scale argument may readily resolve 
this apparent contradiction; the fast time scale of the 
optical measurement should enable one to detect slowly 
fluctuating order parameters as virtually a long-range or- 
der. This situation reminds us of the optical conductivity 
spectra of the well-known CO material a'-NaV20s JL2|. 
The spectra for the CO phase of this system differ very 
little from the spectra of the uniform phase at temper- 
atures above Tool there exists slow fluctuation of the 
charge degrees of freedom even far above the transition 



temperature |l^, Q . Similar situation has also been no- 
ticed for the quasi- ID organic materials @j. 

In the present CuO double-chain systems, if the system 
were exactly at quarter filling, the long-range CO might 
occur with vanishing metallic conductivity. This means 
that if a fine filling-control is made on the CuO double 
chains one may be able to obtain the insulating double 
chains. We therefore would suggest that the experimental 
filling control, e.g., by applying pressure [D| or elemental 
substitution (if it purely changes the filling of the 
CuO double chains), seems to be quite interesting. 



VI. SUMMARY 

We have considered the anomalous charge dynamics 
of the CuO double chains of PrBa 2 Cu40s by focusing 
on the optical conductivity spectra recently observed in 
experiment [Q. 

Our study presented in this paper may be summa- 
rized as follows: (i) We have experimental evidence of 
the presence of strong charge fluctuations, i.e., the sys- 
tem may be in the vicinity of the charge-ordered state. 
Then, one may naturally assume the presence of strong 
intersite Coulomb repulsions in the system, (ii) We have 
therefore adopted an extended Hubbard-type model with 
the intersite Coulomb repulsion near quarter filling and 
treated this in the mean-field approximation. We thus 
have found that the possible charge-ordered state may 
be either of the in-line-type or of the stripe-type charge 
modulations, (iii) We have calculated the optical conduc- 
tivity spectra for the model in the same approximation 
and found that the spectra calculated for the stripe-type 
charge ordering reproduce the experimental features of 
the ~40 meV peak and small Drude weight reasonably 
well but that the spectra for the in-line-type charge or- 
dering does not agree with experiment, (iv) Then, we 
have argued that the stripe-type fluctuations of charge 
carriers should present in the system and this is the ori- 
gin of the anomalous charge response of the system. The 
CuO double chains have no long-range charge ordering 
but the fast time-scale of the optical measurement should 
detect the very slow fluctuation of charges as virtually a 
long-range charge ordered state. 

The CuO double chains in PrBa2Cu40s thus provide 
a good opportunity for studying the anomalous charge 
dynamics in the ID strongly correlated electron systems. 
We hope that further experimental and theoretical stud- 
ies will be made to clarify various aspects of this intrigu- 
ing system. 
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